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ON CROSSED-BEAM MONITORING OF ATMOSPHERIC WINDS AND 
TURBULENCE WITH TWO ORBITING TELESCOPES 

By 

F. R, Krause, S. S. Hu, and A. J. Montgomery 
George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 


This report documents the phase A feasibility study on the MSFC 
inflight experiment No. 23, entitled n Crossed-Beam Arrangements for 
Remote Sensing Missions. n The crossed-beam test arrangement employs 
triangulation of two collimated beams of radiation for the remote sens- 
ing of the space and time variations of thermodynamic state variables 
such as tracer concentrations. The generation, convection and decay of 
trace constituents (aerosols, ozone, and water vapor) will modulate the 
mean radiative power, which is put into the telescope f s field of view by 
an extended radiation background such as scattered sunlight or thermal 
emission. The associated fluctuations of radiative power are then cor- 
related to retrieve local information at the beam intersection from the 
integrated optical signals. The potential of crossed-beam test arrange- 
ment is illustrated by discussing analytically the measurement of 
horizontal wind profiles. It is conceivable to conduct crossed-beam 
missions with two astronomical telescopes (ATM) which are mounted on a 
single spacecraft in such a way that they point towards the earth. 
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TECHNICAL MEMORANDUM X- 53538 


ON CROSSED-BEAM MONITORING OF ATMOSPHERIC WINDS AND 
TURBULENCE WITH TWO ORBITING TELESCOPES 

SUMMARY 


This report documents the phase A feasibility study on the MSFC 
inflight experiment No. 23, entitled "Cross ed-Beam Arrangements for 
Remote Sensing Missions. " The crossed-beam test arrangement employs 
triangulation of two collimated beams of radiation for the remote sens- 
ing of the space and time variations of thermodynamic state variables 
such as tracer concentrations. The generation, convection and decay of 
trace constituents (aerosols, ozone, and water vapor) will modulate the 
mean radiative power, which is put into the telescope* s field of view by 
an extended radiation background such as scattered sunlight or thermal 
emission. The associated fluctuations of radiative power are then cor- 
related to retrieve local information at the beam intersection from the 
integrated optical signals. The potential of crossed-beam test arrange- 
ment is illustrated by discussing analytically the measurement of 
horizontal wind profiles. 


I. INTRODUCTION 


Crossed-beam test arrangements of ground and orbiting detectors 
have been studied to determine whether the remote sensing of horizontal 
wind profiles and turbulence parameters is feasible. Two remote sensing 
devices are used to measure the space and time variations of local emis- 
sion and extinction processes. The fluctuations of the two optical 
signals are combined in a cross-correlation analysis to retrieve the 
signals which are common to both beams out of atmospheric and instrument 
noise in a way known from signal processing and random vibration analysis. 
This combination of remote sensing and cross-correlation analysis has been 
called the M crossed-beam concept." Early versions were proposed inde- 
pendently by M. J. Fisher [1] for jet noise studies and M. Wolf [2] for 
airglow studies. These results showed that correlations from hot wire 
and crossed-beam observations in subsonic jets yield identical results 
for convection speeds [1] and that crossed-beam observations of the air- 
glow altitude are in agreement with rocket observations [2]. Thus, the 
crossed-beam correlation method is theoretically sound, and it has been 
shown experimentally to provide results in agreement with more accepted 
methods . 


In most remote sensing methods, the attenuation or emission of 
electromagnetic radiation is measured along the entire path from the 
source to the detector. To interpret the results of such an observa- 
tion, some assumption is generally required. For example, in atmos- 
pheric studies, one might use the adiabatic or isothermal model 
atmosphere. Assumptions for dynamic studies are made concerning sym- 
metry, expansion properties, and temperature. The important point is 
that, with single path observations, it is not generally possible to 
obtain spatially resolved measurements of the thermodynamic or flow 
properties without invoking such models. 

In crossed-beam experiments, theoretical models need not be invoked. 
The spatial and temporal variation of turbulence and, hopefully, thermo- 
dynamic properties would be measured in sometimes inaccessible regions, 
where other methods are restricted to a mean value over the entire 
optical path or to a direct measurement in a single point. The poten- 
tial applications of crossed-beam experiments, therefore, include almost 
every conceivable experiment in which the basic measurement requires 
monitoring space and time variations of atmospheric species and con- 
taminants. This could include observations of low level and clear air 
turbulence, motion and distribution of air pollution from industrial 
wastes and rocket tests, evaporation over oceans, and ozone motion and 
distribution (which bears some relation to global weather patterns). 

It might also be possible to study the nature and motion of targets of 
opportunity such as storms, hurricanes, meteors, cometary dust, and 
noctilucent clouds. 

The feasibility and potential of crossed-beam experiments in the 
lower atmosphere are illustrated by considering the operation which 
would be necessary to measure a horizontal wind profile up to altitudes 
of 30 km. In chapter III, the theory of space-fixed beam crossings is 
reviewed to indicate how repeated experiments with different beam separa- 
tion distances provide for an accurate statistical description of the 
space and time variations in any random field which produces detectable 
fluctuations of radiative power. Chapter IV gives the ultimate objective 
of in-flight experiments with crossed beams. In chapter V scattering 
and reflection of sunlight, thermal radiation from the atmosphere or the 
earth and non thermal emission in the atmosphere are introduced as extended 
natural sources of radiative power that need not be tracked. Measured 
modulations of laser light indicate that atmospheric fluctuations do pro- 
duce detectable fluctuations, even under a clear sky, which are not 
obscured by receiver noise. The use of these natural sources then pro- 
vides for extremely flexible test arrangements, which could sweep the 
entire altitude profile of horizontal winds with a single fly-by using 
orbital detectors or with a repeated scan, using ground detectors. A 
global sweep requires crossed-beam scanning, the theory of which is 
given in chapter VI. In chapter VII a series of fundamental ground 



experiments is outlined which would ultimately simulate the sweep with 
orbiting detectors. Preliminary specifications are given for both 
ground and orbital detectors. 


II. DEFINITION OF SYMBOLS 
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Symbol 


a 


P 


b 

B 

D 

E 


e 

f 


I 

i = I - 1 


K 


I Si 

I Sx 


k = K - K 
L 

x,y,z 


N 

P 

P 


Q 

R k ( ^ T) i?,A 


Definition 

beam inclination relative to the plane of the sweeping 
beam normal 

inclination of intersection paths relative to beam 
front 

baseline or distance between telescopes 
bandwidths of data acquisition chain in cps 
eddy diffusion coefficient for mass transport 
eddy diffusion coefficient for momentum transport 
charge of one electron 
turbulence frequency 
detected radiative power 

detectable power fluctuation or optical signal 

extinction or emission coefficients (simultaneous 
emission and extinction are to be excluded) 

fluctuations of extinction coefficient 

integral scale of turbulence in x,y,z direction, 
respectively 

number of desired altitude points 
length coordinate along beam trajectory 
path length considered in the averaging procedure 
quantum efficiency [electrons /photon] 

= k (x, t , A) k(x + g, t 4* t, A) two point product mean value 
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Definition 


density of air 

= i (t) i (t + t) correlation of detected fluctuations 

("two beam product mean value") 


T static temperature or integration time 

U,V,W velocity components in space-fixed frame 

W = W + w vertical velocity component 

x horizontal length coordinate normal to both beams 

y horizontal length coordinate parallel to beam front 

z vertical coordinate 

x = (x,y,z) point vector in space-fixed Cartesian frame 

T),£ beam-fixed coordinates relative to the line of minimum 

beam separation 

A optical wave length 

AA wave length interval set by the optical filter 

A/AA spectroscopic resolution 

| minimum beam separation defining wave normal 

— > 

| = (|,t],£) point vector in coordinate system occupied by beam structure 
at time zero 


R(0,0W , 


space time correlation coefficient 


time lag 


Superscripts 


vector 


ensemble, space and/or time average 
time average for moving observer 
length average over intersection path 
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Subs cr ipts 

v z,y or 1,2 line of sight in z,y direction 

* moving frame autocorrelation 

C group velocity component (convection of disturbances 

normal to beam front) 

A,B points of minimum beam separation 

III. REVIEW OF THE CROSSED-BEAM METHOD 


1 . Optical Signals from A Single Beam 

In a crossed-beam arrangement the local flow region of interest 
may be chosen by triangulation between two collimated beams and extends 
between the two points x and x* + (|,0,0) of minimum beam separation. A 
typical test arrangement for jet noise studies is shown in figure 1. The 
beam diameters (viewing angle) and orientation (azimuth and elevation) 
are set by a remote telescope. The collected radiative power is then 
filtered by a monochromator, which is set at the wave length A and the 
resolution AX as illustrated in figure 2. The remaining signal, 1^, is 
then monitored by a photodetector with a fast response such that the time 
history is not lost. 

The remote sensing device will see any local emission or extinc- 
tion process inside the beam. We choose to express the contributions to 
the radiative power AI, which originate in a beam element, Ax, by a 
generalized spectral extinction coefficient 

K(?,t,A)-ig. (1) 

The word "generalized" indicates that emission processes are included. 

As long as either emission or extinction is locally dominant, the equa- 
tions of radiative transfer show that the total signal, I, may be expres- 
sed in terms of a line integral over the extinction coefficient. 

detector 

I (t, A) = I exp (- / K(x,t,A) dx) . (2) 

d source \ J / 

source 
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The fluctuation of this signal, 


1 - x d - ij’ < 3 > 

can be calculated by subtracting a mean value I, which may be established 
by averaging over time and/or over repeated experiments. In case of 
ideal detectors, this fluctuation is related to the fluctuations of the 
extinction coefficient 


k = K - K 


( 4 ) 


through the line integral of equation (2). 
detector 

(K + k) dx 

source 

. = I e 

d source 



= (I s + i s ) exp 


detector 



source 



(5) 


It is now assumed that the integrated k fluctuations of the medium can 
be treated as small perturbations. 


detector 

k(x y t y 7\) dx | «1. (6) 

source 



In this case the last equation can be rearranged to give the desired 
fluctuation: 


i = I, - I, = — i - l, / 

d d - source d J 

s source 


detector 

f k(x y t,7\) dx. 


(7) 
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The first term gives the contribution of the light source noise, and the 
second accounts for the fluctuations of the medium. Detector noise has 
to be added for nonideal detection devices. 

The small perturbation assumption is not very restrictive. The 
transmiss ion 


J K dx, 

as well as the local values of the fluctuation k could be large. Only 
the fluctuation integral over several statistically independent parts 
has to be small enough to permit linearization. 

2. Retrieval of Local Information 

The way the beam correlation works to retrieve the local informa- 
tion at the region of minimum beam separation becomes clear by indtroduc- 
ing the concept of a "correlation volume." This volume is defined in 
turbulence theory by an imaginary experiment with two-point observers, 
as illustrated in figure 3. We consider a special arrangement where the 
beams are aligned parallel to the y and x axis and cross in the space- 
fixed position x. The first imaginary observer shall be attached to this 
position. The^ second observer monitors the k fluctuations at the variable 
position 5? + §, and the imaginary experiment now calls for multiplying and 
averaging the observed fluctuations. This gives a "two-point product mean 
value," 




* > r ~ i ? 

,t,A) = k(x,t,A) ic(x + t + t, 


A; : 


r q 


) 


which is the basis of all analytical descriptions of random fields [5]. 
The random nature of a fluctuation field such as k (5?, t,A) is recognized 
by the fact that the two-point product mean value drops to zero over 
finite distances |1*j . This leads to the definition of the correla- 
tion volume. It is that volume surrounding a certain point x , within 
which the two-point product mean value R^. remains finite. The average 
radius of this volume in the x-direction is used to define the integral 
length scale: 


L 

x 


.L 

2 



-00 


k(x + fc, v, z, t) k(x,t) 
k 2 (?,t) 


d|- 


(9) 
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Similar expressions hold for the extensions 2Ly and 2L Z of the correla- 
tion volume in the y and z directions. In wind tunnels typical L values 
are between 1 mm and 10 cm. In the atmosphere they vary between 10 m 
and 1 km. 


A crossed-beam experiment resembles the above imaginary experi- 
ments by performing exactly the same calculations. Only the local fluc- 
tuations of k are replaced with the fluctuations i of the detected 
signal, I. This leads to the M two-beam product mean value 11 


R(S,t,A) - i y (t, A) i 2 (t + t, A) . (10) 

The relation between two-beam product and two-point product mean values 
is now derived by splitting the detectable optical fluctuations iy and 
i 2 into three different components, as shown in figure 3. 


i = i + i + i 

y yi y 2 y3 


The first component describes the contribution of the correlation volume 


1 . = 

yi 


k (x + i, y + tj, z, t + T, A) dr). (11) 


correlation volume 


The second contribution accounts for the medium fluctuations outside 
the correlation volume; that is, for the flow noise. 


detector 


ys 


k dr] - i 

1 yi 


( 12 ) 


source 


The third contribution summarizes light source noise and detector noise, 

I „ 


i = i h detector noise. 

y3 s - 

s 


( 13 ) 


c 
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The ins tantaneous product of the two optical signals may thus be xcrritten 
as the sum of nine other products, 


ii=i i + i i 4- . . . i i . 
y z yi zi y 2 Z 2 y3 z3 


(14) 


The last eight of these products have one noise component as a factor. 
Their instantaneous value will thus assume positive and negative values. 

If we average these products, their sum will decrease to zero with at 
least the inverse square root of integration time and/or number of experi- 
ments. The first product contains the fluctuations inside the correlation 
volume which are common to both lines of sight. The instantaneous product 
of these fluctuations is always positive, and its average value will thus 
approach a constant value with increasing integration time and/or number 
of experiments. This product will outweigh all others if we add long 
enough. In the limit of infinitely long summation, one should be able 
to retrieve almost infinitesimal contributions iy 1 and i Zl . 


R = i (t) i (t + t) = i i = A)-> (15) 

z v y 7 yi z i * 'x 

Accordingly, the cross-correlation has pulled the correlated signals from 
the correlation volume out of the flow noise (i y2 , i 22 ) anc * instrument 
noise (iy 3 ,i z3 ). Substituting the line integrals for iy 1 and i 2l , equa- 
tion (11), inverting the order of integration and averaging, gives the 
following alternate form of the last equation: 


R(I,t,A)^ = I y i z J' k (x,y, z - (^t.A^kCx + y + r\, z, 

correlation 

volume 

The crossed-beam method can do better than resolving only the 
tion volume. This follows from a more detailed consideration 
optical integration. The above product of two line integrals 
transformed to an area integral, 


R(!,t,A) ? = ll z JJ R k (‘?,T,A) i? d T1 d£, 

correlation 
"disc" 

if the fluctuations are locally homogeneous along one beam. The transform- 
ation is identical with the condition of statistical homogeneity, and its 
application leads to an integrand, where only the second point is variable. 

9 


t + t,A 2 ) dT)d£. 
(16) 

correla- 
of the 
can be 


(17) 


k(x,y,z - k (x + %, y + t \, z , t + t, A 2 ) = 


= k(x,y,z, t, Ajl) k(x + y + r\, z + t + t, A 2 ) . 


( 18 ) 


The area integration is illustrated in figure 4. The integrand is equal 
to the two-point product mean value between two observers. The first 
observer is fixed to the upstream point x* of minimum beam separation. 

The second observer may be anywhere in a plane which is parallel to both 
beams and contains the second point (x + g, y,z) of minimum beam separa- 
tion. This plane will be called the beam front. The two-beam product 
mean value is thus identical with an optical integration over the beam 
front. Since along the beam normal no integration takes place, it is 
possible to resolve two-point product mean values ins ide the correlation 
volume by repeating the experiment for several beam separations, |. 

The following analytical considerations, as well as M. J. Fisher 1 s 
experiments [6], lead us to the conclusion that the area integrations are 
weighing the points close to the beam normal such that the normalized two- 
beam product mean value 0 closely approximates a two-point product mean 
value between the two points x and x + (|,0,0) of minimum beam separation: 


R(i,T,A)^ tyi.o.o, T, A) i y U * 0) 
0(£,t,A)£= r( 0,0,A) ? " 1^(0, 0,0,0, A) i (| = 0) 


(19) 


The errors of this approximation are small in the two extreme cases of 
very small and very large lateral scales Ly and L z . In case of very small 
lateral scales the weighting of points close to the wave normal is self- 
explanatory. In case of a large lateral scale, the fluctuations will be 
provided by an almost periodic wave of infinite extent, that is, in the 
limit by a plane acoustic wave. This follows from the well known property 
of Fourier transform pairs, which states that if the correlation is broad 
(large turbulence scale) then the spectrum (extension in wave numbers) 
will be narrow. However, in case of acoustical waves, a separation-of- 
variables assumption holds for separations in time and along the wave 
normal . 


k(x,y,z,t) k(x+£, y +r), z + f;, 


t+x) = 0 (|, t) k(x,y,z,t)k(x,y + r\ t z+£,t). 


( 20 ) 
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Substituting the separation of variables into both (18) and 
(19), we find that these conditions are met identically. 

For intermediately sized correlation volumes, the two-beam 
approximation of two-point product mean values, equa. (19), is not 
always valid. For locally isotropic turbulence which is closely 
approximated by atmospheric turbulence at the ground, the derivations 
given in the appendix show that the two-point product mean value should 
be directly proportional to the second derivative of the two-beam pro- 
duct mean value with respect to the separation distance: 


isotropic 


Rk<6; T) = 0; i = °; 


x) i?,A' 875^ R(5 ’ t) ?,A 


( 21 ) 


For small separation distances, the function R(f=, . ..) can be approxi- 
mated by a cosine curve, since R(|, ...) is an even function of g 
because of the isotropy. However, the second derivative of a cosine 
curve is directly proportional to the cosine curve itself and the pro- 
portionality factor drops out when normalizing as required by equation 
(19). For small separation distances and locally isotropic fields, 
equation (19) should thus be valid for correlation volumes of any size. 
For large separations, and anisotropic, inhomogeneous turbulence errors 
are to be expected. Fortunately, the experiments discussed in section 
III, 4 indicate that the approximation of point measurements may still be 
good in many anisotropic and inhomogeneous fields of practical interest 
such as jet streams. 

3. Measurable Turbulence Parameters 


The crossed-beam method approximates two-point product mean 
values, and the two-beam product mean values can therefore be used to 
approximate all turbulence parameters which are commonly derived from 
two-point product mean values [8]. 

The local integral scale of turbulence is defined by averaging 
the area under the space correlation, R(g, t = 0), and its graphical 
evaluation is illustrated on figure 5. This scale L x is a measure for 
the extent of correlation volume along the beam normal. For mean value 
estimates as well as error discussions, we would also like to know the 
lateral integral scales L y and L z . These could be measured by rotating 
the lines of sight so that they can be separated parallel to the y and 
z directions. 


V 
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Plotting the beam correlation coefficient against the time 
separation instead of the space separation (figure 6) allows group 
velocities to be read. For a given space separation, the correlation 
will reach a maximum at the time that the turbulent eddies need to 
travel from the upstream to the downstream beam. The time lag, t*,, 
which is indicated by the maximum is therefore a direct measure of group 
velocities. A crossed-beam experiment is thus capable of measuring the 
component of the group velocity, which is normal to the beam front. In 
atmospheric applications, these components will be horizontal, since the 
beams are mostly parallel to a vertical plane. 

A group velocity represents an average velocity of a disturb- 
ance [11]. It will in general differ from a mass average velocity. In 
spite of this fact, the name "bulk convection speed" is often used in 
turbulence investigations [8]. Later sections of this report will use 
the assumption that the group velocity of a crossed-beam experiment 
closely approximates a mass average velocity, i.e., a horizontal wind. 

The two velocities will be equal, if positive and negative amplitudes of 
the correlated signals iy 1? i z i are equally likely, that is, if the first 
order probability density is symmetrical. Extreme cases of unsymmetr ical 
or skewed first order probability densities occur in shear layers. For- 
tunately, the experiments shown in figure 8 indicate that the difference 
between group velocities and mass average velocities is not too great 
even in such extreme conditions. We therefore believe, that a crossed- 
beam system gives horizontal winds with sufficient accuracy to be of 
practical interest. In case of excitation problems, such as a rocket 
going through clear air turbulence, one needs to know the group veloc- 
ities and not the mass average velocities. 

The envelope to the individual temporal correlations (figure 6) 
is called the moving axis autocorrelation, since the associated two- 
point product mean values are those which would be seen by two observers 
moving with the group velocity [9]. The study of the envelope enables 
us thus to study the generation and decay of targets of opportunity such 
as shock waves and weather fronts independent from the convection proces- 
ses. The decay of this correlation curve of 1/e of its original value is 
used to define the average "eddy" lifetime. The Fourier transform of this 
envelope would indicate the dominant frequencies which are felt by the 
moving observer, traveling with the target. 

An observer traveling with the bulk convection speed will record 
almost the same transport phenomena [8] as the observer who travels along 
in a single fluid particle. The moving axis autocorrelation or envelope, 





r < o ’ o >x-,a y* * °> 


( 22 ) 
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is the closest approximation of the Lagrangian autocorrelation functions 
which is accessible to direct or remote sensors. The moving axis time 
scale, 


00 

= f 0 *( t m ) dr M , (23) 

o 

therefore closely approximates the Lagrangian time scale. This allows 
an important contribution to the experimental analysis of turbulent 
transport phenomena. The observed extinction coefficient fluctuations, 
k, are very often directly proportional to the fluctuations p t in the 

tracer concentration which in turn indicate a turbulent mass flux, pj- r w 
across horizontal planes. From the point of view of "turbulent diffusion," 
this mass flux should be directly proportional to the vertical gradient 
of the mean tracer concentration: 


P. w 
tr 


dp. 


= D 


tr 


space- fixed 


(24) 


However, describing the diffusion process through the random walk of a 
single particle, we find that the eddy diffusion coefficient, D, is 
equal to the product between the mean square value of the vertical 
velocity fluctuations and the Lagrangian time scale [10]. 


D 



moving axis 


L 


T* 


(25) 


Combining the last two equations leads to a direct relation between the 
moving axis time scale, L T , and the density gradient: 



n w 
M tr 


space- fixed 


w 


moving axis 



(26) 


The concept of turbulent diffusion can be applied to momentum diffusion 
instead of mass diffusion. In this case, the turbulent flux is given by 
the apparent turbulent shear stress, p uw , and the driving mean velocity 
gradient du/5z: 


— p Bu 

P uw = PE ^ . 


(27) 
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If the eddy diffusivity, E, is approximately equal to the mass 
diffusion coefficient, D (turbulent Lewis number E/D = 1), equation (25) 
can still be used. Its substitution leads to an expression for the 
shear: 


du = 

bz 



space-fixed 




moving axis 



(28) 


Experimental evidence suggests that one-point correlation coefficients 
between velocity components like the one shown on the right-hand side 
are rather insensitive to position on the turbulent field [8]. We might 
therefore expect that the shear is almost inversely proportional to the 
moving axis time scale. In subsonic jet shear layers, this has been 
demonstrated already through the hot-wire experiments of Davies and 
Fisher [7], 


4. Experimental Status 

The first crossed-beam experiments were designed to verify the 
approximation of two-point product mean values in very general fields of 
turbulence. The "imaginary" point probes were approximated by hot-wire 
probes. Although a small hot wire gives velocity fluctuations rather 
than extinction coefficient fluctuations, a comparison between hot wires 
and crossed-beam correlation coefficients should give identical results 
as long as we are interested only in the dynamical description of turbu- 
lence in terms of scales, the group velocities, and the spectra. The 
initial portion of an axis ymme trie, free shear layer was chosen for the 
experiments since the turbulence is anisotropic and inhomogeneous and 
since group velocities, turbulence scales, and eddy lifetimes are well 
documented in the literature. 

All measurements were taken in a subsonic (M = 0.2) air jet 
exhausting through a one-inch diameter nozzle into the atmosphere. 

Light extinction was achieved by spraying a small amount of liquid 
nitrogen into the settling chamber. This produces small water droplets 
in the exhaust flow which attenuate the crossed beams by scattering. 
Later, the nitrogen injection was replaced with water spray, which 
offers a better control of the tracer concentration. 
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A typical sample of optically measured two-beam product mean 
values is shown in figure 7. The beam intersection was adjusted to 
the center of the free shear layer three exit diameters downstream of 
the nozzle lip. The turbulent convection process is clearly indicated 
by the displaced maxima for finite beam separation. The envelope to 
these maxima falls to 1/e of its initial value for a time lag of 
880 microseconds. This is in excellent agreement with the hot-wire 
data from which we estimate that the value of the eddy lifetime for 
our experimental conditions should be 920 microseconds. 

Figure 8 shows the same experiment under different conditions. 

The initial intersection point of the beams was six diameters from the 
jet exit, (i.e., at x/D = 6.0) instead of at x/D = 3.0). Instead of 
using liquid nitrogen to generate a water droplet fog tracer, a small 
nozzle was mounted in the jet settling tank which sprayed water droplets 
into the flow to form the tracer fog. The curves appear to fall onto 
or close to an envelope as required. The moving axis autocorrelation 
curve has been obtained from the hot-wire data of Davies, et al., [7] 
and is shown as the dotted line superimposed on the measured cross cor- 
relation curve. 

The spatial resolution of the crossed-beam method seems very 
encouraging, considering that the envelope from the hot-wire data 
exactly matches the individual curves from the crossed-beam me thod . 

For the moving frame, the optically measured autocorrelation agreed 
with the hot-wire curves, regardless of whether nitrogen or water spray 
was used as a tracer. 

The optically measured results in figure 8 were only 0.2 inch 
apart. This is much smaller than the geometrical shear layer thickness 
.8 in) which is indicated by the hot-wire results. Thus, space-time 
correlations have been resolved inside a correlated volume by moving 
the intersection point. 

The spatial resolution obtainable with crossed-beam test 
arrangements becomes more apparent by comparing optically measured 
group velocities with hot-wire group velocities. Our main results are 
summarized in figure 9. In the central and low speed portions of the 
shear layer (q >0), the convection velocity determinations agree within 
experimental accuracy. For q < 0, however, the crossed-beam results 
tend to follow the mean velocity profile (U/U Q ) rather than attain the 
constant maximum value observed for hot-wire data. The experiments of 
Davies [12] indicate that hot wires at the high velocity edge of a shear 
layer react to sound radiation from the shear layer center as well as to 
the convection of upstream turbulence. A low value of the hot-wire group 
velocity is therefore to be expected because of the convection of the sound 
sources which travel in the middle of the shear layer with lower velocity 
than the edge particles. Instruments other than the hot wire should give 
higher group velocities, as indicated by the crossed-beam method. 
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The optical group velocity measurements used tracers. Since 
then, we have succeeded in measuring convection speeds in supersonic 
jets without tracers. These tests used the density fluctuation of air 
to produce changes of radiative power in the ultraviolet absorption 
bands of oxygen (Schumann-Runge bands). Since the method worked on 
both scattering and gaseous absorption, two-phase flow application 
became conceivable such as condensation and evaporation of water vapor, 
or the spreading of rocket exhaust products, and unburnt fuels. 

Ultraviolet absorption measurements with collimated beams 
proved to be extremely difficult because of the comparatively low levels 
of the mean radiative power. The commercially available Hannovia lamps 
and McPherson lamps are much less bright than sources emitting visible 
light. Also, the small beam diameter (2mm) severely reduces the solid 
angle that can be used to collect radiation from the source. There- 
fore, we had unusual low power levels, that is, high shot-noise levels. 
The usual method of reducing shot noise, using electronic filtering, 
could be employed only to a limited extent since the bandwidth of the 
filter must be such that the complete turbulent spectrum is accepted. 
Thus, the shot noise was one order of magnitude higher than the 
RMS level of the integrated optical signal. In spite of these limita- 
tions, nine out of twelve runs showed the correct convection speed. 

Later tests used a focused beam system rather than a collimated beam 
system [13]. This increased the signal to noise ratio by a factor of 
200 . 


IV. CONCEIVABLE IN-FLIGHT EXPERIMENTS 


We propose that crossed-beam test arrangements are used in remote 
sensing missions to monitor horizontal winds, humidity and thermodynamic 
properties of the earth’s atmosphere [14]. Two typical test arrange- 
ments are shown in figure 10. One illustrates the local use of ground 
detectors and the other one a global coverage with orbiting detectors. 

Present meteorological sensors are mounted on meteorological towers, 
balloons, aircraft, and rockets. These monitoring systems have the 
following inherent limitations: 

(1) Single path experiments, 

(2) Restricted altitude ranges, 

(3) No control over horizontal traverse, 

(4) Interference with the phenomena of interest, 

(5) Large-scale averaging. 


I 
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Because of these limitations, some meteorological problems partially 
defy the use of existing instrumentation [15]. Meteorologists have 
already turned to remote sensing devices to overcome part of these 
instrumentation problems. It is hoped that remote sensing devices 
accomplish the following three objectives [16], 

(1) Establish global observation of the entire earth's 
atmosphere on a regular basis. 

(2) Provide observations for direct use in mathematical 
prediction models. 

(3) Follow the full life cycle of weather systems such as 
hurricanes . 

In meeting these objectives, one could overcome the two main problems 
which at present plague numerical weather prediction [16]; (1) insuf- 

ficient grid points and (2) unknown effect of small scale disturbances 
on the nonlinear behavior of larger scale motions. 

At present a grid of observation stations is available only in 
North America, Europe, and the northern part of Asia. Except for 
Australia and very few observation stations in Africa, South America, 
and Antarctica, the Southern Hemisphere is void of data. Attempts are 
now being made to measure winds with long-lived constant-pressure 
balloons, tracked either from the ground or by satellite. The spatial 
positions of these balloons depend on the circulation pattern and not 
on the desires of the meteorologist/ Any improvement in providing 
meteorological information at discrete levels in the atmosphere (« 50 
to 100 mb pressure level), at discrete time intervals (« 4 times a day), 
and with discrete spacing 200 miles) would be greatly welcomed and 
would be considered a big leap forward in improving forecasting condi- 
tions. The remainder of this paper indicates how crossed-beam test 
arrangements may be used in remote sensing missions to supply the mis- 
sing information at the desired location. 

The second problem was to estimate the effect of small scale disturb- 
ances, such as convective cloud patterns, small vortices, turbulence, 
etc., on the nonlinear behavior of larger scale motions. Attempts are 
now being made to estimate such small scale systems from a stationary 
satellite which takes cloud pictures over the Pacific Ocean. However, 
commensurate information on winds, temperature, etc., is still lacking. 
This means that, in addition to the large scale information mentioned 
previously, small scale measurements should be available for certain time 
intervals and certain critical regions, which may be specified on short 
notice by weather prediction centers. Here again, a crossed-beam moni- 
toring system might be the answer. On the following pages we will 
indicate that scanning crossed-beam systems give turbulence length 
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scales and wave number components as well as wind information. These 
results are necessary to study the transfer of kinetic energy between 
atmospheric motions, which are characterized in terms of macroscales, 
mesoscales and microscales. Single-point probings such as balloon-borne 
beam instruments cannot provide this information. 

One other important application would be an early warning system 
for clear-air turbulence. A crossed-beam arrangement with optical detec- 
tion devices should be able to detect turbulence in atmospheric regions 
which are "clear 3 * * * * * * * 11 to radar receivers. A space-fixed system might be 
used to monitor clear-air turbulence above traffic centers, whereas an 
orbiting system could be used to sweep desired flight paths. 


V. EXPECTED OPTICAL SIGNALS 


The feasibility and potential of crossed-beam experiments in the 
lower atmosphere are illustrated in this paper by considering the opera- 
tion which would be necessary to measure a horizontal wind profile up to 
altitudes of 30 km. Such crossed-beam arrangements, figure 10, would 
work whenever detectable space and time variations of emission or extinc- 
tion processes are common to both beams. This requires local changes in 
radiative power of sufficient intensity that the change may still be 
resolved after being transmitted to the distant spectrophotometer. 
Furthermore, the radiation source should be extended rather than a 
point source to avoid fluctuations caused by refractive index variations, 
such as the scintillation of star images. We have therefore tried to 
answer the following questions: 

(1) What is the mean radiative power of suitable extended 
radiation sources? 

(2) What modulations of this radiative power can one expect 
due to atmospheric inhomogeneities which move broadside 
relative to single beam? 

(3) Are these meteorological modulations above the detector 

noise levels? 

1. Mean Radiance of Extended Natural Sources 


Consider a photometer coupled to an astronomical telescope. 

The detectable radiative power is contributed by all radiative phenomena 

inside the collimated field of view of the telescope (see figure 11). 

In atmospheric applications, the detector is placed on either the ground 

or the orbiting spacecraft as shown in figure 10. Because of light 

absorption in the atmosphere, remote sensing of meteorological phenomena 
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inside the atmosphere requires work in regions of the electromagnetic 
spectrum, where the atmosphere is at least partially transparent. 

These "optical windows" may be found from typical transmission spectra 
like the one shown in figure 13. The bottom part of this figure shows 
the spectral region where suitable windows exist. 

Remote sensing within the "optical windows" requires the exist- 
ence of either man-made or natural radiation sources. The use of man- 
made sources would severely limit the application of the technique, and 
point sources such as stars must also be discarded because of scintil- 
lation effects. Natural, extended background sources which require no 
tracking and suffer no scintillation are therefore ideal sources for 
use in the crossed-beam technique. 


The following sources are conceivable: 


(1) Scattered solar radiation 

(2) Emission from the atmosphere 


ground detector 


(3) Solar radiation reflected and 
scattered by the earth* s sur- 
face 

orbital detector 

(4) Thermal radiation from the 
earth and the atmosphere 


All these sources constitute an extended background of radiation. We 
will discuss the characteristics of these sources that are potentially 
useful as part of a weather monitoring system. The main results are 
summarized in table I. 


Ground-based detectors (figure 10) can use scattered sunlight 
and infrared emission from the atmosphere as radiation sources to 
monitor atmospheric phenomena. 

Scattered Sunlight . Solar radiation is both absorbed and 
scattered in passing through the earth* s atmosphere to an extent which 
varies with the wave length. The earth's surface receives radiation 
both in the form of a parallel beam of rays from the sun and scattered 
radiation incident from all points of the sky. The total flux of 
scattered radiation reaching an object depends on the sun angle, the 
cloud cover, the terrain reflectance, and the concentration of tracers 
in the atmosphere. The flux of the scattered radiation is weakly 
dependent upon earth latitude and more strongly dependent upon sun 
elevation angle, as shown in figure 12. 


19 


Table 1. Estimation of Radiance Available to 
Ground and Orbital Detectors 


Detector 

Location 

Ground 

Ground 

Orbit 

Orbit 

Extended 

Source 

Scattered 

Sunlight 

Ozone 

Emission 

Reflected and 
Scattered Sun 
Light 

Thermal Emis- 
sion of Earth 

Monochromator 
Bandpass AA 

0.38}i to 0.72jj, 

9p. to 10|i 

0. 38jj, to 0. 7 2jj. 

8p. to 12)j. 

Estimated 

Mean 

Radiance I 

4.5 • 10" 3 
watts/ cm 2 s ter 

7.8 • 10" 5 
watts/ cm 2 ster 

4.2 • 10“ 3 
watts/ cm 2 ster 

4 • 10 -3 
wa t ts / cm 2 s ter 

Estimated 
Modulation 
in Spectral 
Radiance i^ 

4.5 • 10“ s 
watts/ cm 2 s ter 

7.8 • 10~ 8 
watts/ cm 2 ster 

4.2 • 10" 6 
watts/cm 2 ster 

4 • 10“ 6 
watts/ cm 2 ster 

Modulation 
Bandwidth B 

300 cps 

300 cps 

1000 cps 

1000 cps 


Influenced by the various parameters mentioned above, scattered 
solar radiation fluxes have been measured ranging from 0.35 to 2.80 x 10* 2 
watts/cm 2 [17]. The same reference reports individual cases where the 
values of the flux of the scattered radiation exceeds 7 x 10~ 2 watts/cm 2 . 

A conservative figure for the flux from a cloudless sky with a solar alti- 
tude of 45° would be about 1.4 x 10~ 2 watts/cm 2 . 

In crossed-beam experiments we are not interested in the total 
flux received from the entire sky (hemisphere = 2 it steradian), but only 
in the sky radiance which is available to the small angle of view £a, of 
a telescope [4 it sin 2 (Aa/4) steradian]. 

Assuming that Lambert’s law gives a reasonable approximation of 
the angular distribution of scattered sunlight, the radiance measured by 
a vertically pointing telescope becomes 


1.4 • 10’“ 2 /tc = 4.5 • 1CT 3 watts/cm 2 steradian. 
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Of prime importance in our study of the scattered radiation is 
the radiance available in an emission bandwidth, AA. The spectral 
radiance of scattered radiation is shown in figure 13 for a cloudless 
sky with the sun at a zenith angle of 45° for an observer located on 
the earth 1 s surface and at an altitude of 30 km [15], Taking account 
of detector response, a useful bandwidth of emission for scattered 
solar radiation is from 0.38 to 0 . 7 2 pi with a maximum intensity between 
0.4 and 0.5p,. Radiation detected at wave lengths greater than about 
2 - 3 pi is not due to scattering of solar radiation, but to thermal 
emission in the atmosphere. 

The flux of scattered sunlight is, in first approximation, 
roughly proportional to the density of the scattering molecules above 
the detector. At 30 km the atmospheric density is 10~ 2 , that of sea 
level. Hence, we may expect the local source strength of scattered 
solar radiation at these altitudes to be of the order of 10 2 times 
smaller than at sea level. This is found to be so, as reported in [15] 
and shown in figure 13. 

Ozone Emission . Besides scattered sunlight, the ground 
detector could also detect the thermal emission from the atmosphere. 

For crossed-beam experiments we need an atmospheric window to n see" 
through the troposphere. However, there has to be some emission, 
absorption, or scattering of radiation present for application of the 
crossed-beam technique. In this study therefore, the feasibility of 
application of the technique using the 9.6 thermal emission band of 
ozone as the source of radiation was considered. Since this radiation 
band falls within the 8 p - 12p atmospheric window, measurements through- 
out the troposphere are possible. This band is clearly indicated by the 
"dip" in the atmospheric transmission (figure 13). For ground dectectors 
this band is important because: 

(1) It is in a spectral region where it is the only 
active tracer; there exists no water vapor or 
strong C0 2 bands to interfere. 

(2) Contributions to the signal at 9.6p, received by a 
ground detector from altitude ranges 0-10 km, 10-24 km, 
and above 24 km are at the same order, and thus ozone 
emission may be used in crossed-beam experiments. 

The following table shows the spectral radiance of ozone emission in the 
9.6 fj. band as a function of zenith angle, with the background continuum 
subtracted [18]. 

The atmosphere was arbitrarily divided into three layers, and 
the contribution from each is given. 
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Table 2 


9. 6p Ozone Emission Intensity (watts cm” 2 ster" 1 p" 1 ) 


Zenith Angle 

Ground to Tropopause 
(0-10 km) 

Tropopause to 30 MB 
(10-24 km) 

Above 30 MB 
(above 24 km) 

0 

1.6 x 10 -5 

2.7 x 10" 5 

3.5 x 10" 5 

60° 

3.0 

3.8 

4.0 

75 °32' 

5.2 

5.2 

3.5 

80°24' 

7.2 

5.9 

2.8 

82°48' 

8.7 

6.0 

2.6 


The first line of Table 2 provides an estimate of the total 
radiance which was measured by a ground detector looking vertically 
upward through the atmosphere at a wave length of 9.6 microns. If this 
ozone emission band is taken to be one micron wide, i.e # , AA = lp, then 
the total detected intensity will be 


7.8 x 10" 5 watts cm" 2 ster*" 1 p” 1 . 


More recent measurements by Bell et al. (1960) reproduced in 
the "Handbook of Geophysics and Space Environments" AFCRL, give higher 
figures for the spectral radiance of the ozone emission. For the 
zenith direction, the radiance at the middle of the ozone band was 
measured as 2.5 x 10" 4 watts cm" 2 ster" 1 p" 1 , and for a zenigh angle 
of 60°, 4.5 x 10” 4 watts cirT 1 ster" 2 p" 1 was recorded. However, we 
have chosen to use the more conservative number above. 

Orbital detectors may also detect scattered sunlight and 
ozone emission. Additional extended radiation sources are provided by 
reflected sunlight and by the earth's thermal emission. The following 
discussions are restricted to the same atmospheric windows chosen for 
discussion of ground detectors, although application of the crossed- 
beam technique to other wave lengths are in no way ruled out. 
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Scattered and Reflected Sunlight . The radiation of wave length 
A received by a detector looking toward the earth from an orbiting satel- 
lite depends on a large number of factors: 

(1) Intensity of sun's radiation incident at edge of 
atmosphere at wave length A. 

(2) Transmission of atmosphere. 

(3) Zenith angle. 

(4) Contribution by scattering from molecules and aerosols. 

(5) Diffuse reflectance of terrain. 

(6) Emission of earth and atmosphere. 

(7) Presence or absence of clouds. 

(8) Spectral bandwidth of detecting system. 

(9) Collector optics. 

(10) Detector sensitivity. 


The intensity of the sun's radiation at the edge of the atmos- 
phere and at ground level and diffuse reflectances of different types of 
earth terrain are given in table 3. 


Table 3 

Spectral Radiance of Sun at Normal Incidence 


Outside Atmosphere 


Wave length (p) Intens ity watts/cm 2 / p 


0.4 

0.5 

0.6 

0.7 


1.74 x 10" 1 
1.96 
1.8 
1.5 


At Sea Level-Zenith Angle 60° 


Wave length (p) Intensity watts /cm 2 / p 


0.4 

0.5 

0.6 

0.7 


0.47 x lO -1 
1.2 
1.17 
1.11 
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Table 3 (Continued) 

Typical Values of Terrain Reflectance (0.4-0. 7u) 


Terrain 


Reflectance 


Snow 0.8 

Water 0.05 

Forests 0.03 

Open grassland 0.06 

Dark brown soil 0.11 


On the basis of these figures, the radiance available to orbital 
detectors from the sunlit portions of the earth's surface may be estimated. 
Taking a terrain reflectance of 0.25, the contribution by reflection from 
the surface will be of the order of 


1.8 x 10 


-l 




x (0.7) 2 x 0.25 = 7 x 10~ 3 watts cm" 2 ster" 1 micron" 1 


= 2.11 x 10" 3 watts cm" 2 ster" 1 


for the spectral bandpass 0.4 — 0.7 p.. A mean atmospheric transmission 
of 0.7 is assumed here, corresponding to a zenith angle of the order of 
30°. 


In addition, there will be a contribution due to scattering by 
molecules and aerosols in the atmosphere which will be strongly dependent 
upon the direction of observation, direction of the sun, and also the 
ground reflectivity. Calculations have been made for a plane parallel 
Rayleigh atmosphere by Coulson, Dave and Sekera [19] for a large variety 
of conditions. Assuming a zenith angle of 30°, a ground reflectivity 
of 0.25 then from tables the radiance seen by a detector looking verti- 
cally downward through the atmosphere in the wave length band 0.4 - 0.7 
microns will be 4.2 x 10" 3 watts cm" 2 ster" 1 . This is a lower limit and 
was used in table 1. The presence of aerosols will increase the scatter- 
ing contribution and decrease the contribution of reflection of the earth's 
surface. An upper limit is obtained by assuming a diffuse reflectance of 
unity, which gives 1.7 x 10~ 2 watts cm" 2 ster" 1 for the same spectral 
bandwidth. 
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Thermal Emission from the Earth . Besides reflected and scattered 
sunlight, the orbital detector could use the thermal emission from the 
earth as an extended radiation background. Most substances on the earth's 
surface have nearly unit infrared emissivity. In the relatively trans- 
parent infrared regions of the atmosphere and on clear days, this radia- 
tion is nearly free to escape unaltered into space. The main variations 
observed are from equator to pole as the earth's surface temperature 
changes. The radiance of the sun side of the earth observed outside the 
earth's atmosphere by a satellite is shown in figure 14. In the trans- 
parent regions, the radiance corresponds to a black body at the surface 
temperature of the earth, while in the strongly absorbing regions, it 
corresponds to the black-body radiance at stratospheric temperatures. 
Between 8 and 12 microns, a satellite-borne detector will therefore see 
the same radiance both on the night and on the day side of the earth. 

Since the 8 to 12 micron band of emission is near the peak of the 
black-body radiation, the thermal emission, which is picked up by an 
orbiting telescope, can be calculated from the black-body curve. The 
associated radiance is 

4 • 10“ 3 watts /cm 2 ster. 


2. Meteorological Power Modulations 

Meteorological modulations of the mean radiative power are pro- 
duced by the generation, convection and decay of emitting, scattering 
and absorbing particles which follow the atmospheric motions. The follow- 
ing discussion is restricted to scattering. The only direct measurements 
of such extinction modulations currently known to us have been made in 
conjunction with optical communication studies. 

Meteorological modulations of a He-Ne laser have been measured 
over optical paths up to 6.8 km long (mountain range) [20]. The laser 
output of 0.8 milliwatt was collected by a 24- inch reflecting telescope. 

By using the full aperture of the telescope, no "spill-over" of the laser 
light occurred and the attenuation measured by the detector was due to 
atmospheric attenuation and not to refractive index variations which 
would cause the bending of the light beam away from the sensor. Figures 
15a and 15b give the power spectrum of the intensity fluctuations. The 
integration, which is outlined in figure 15a, indicates that the ratio of 
the root-mean-square value to the mean (d.c.) level of the fluctuating 
signals are on the order of one percent. This one-percent fluctuation 
is taken as representative of the intensity fluctuations occurring over 
large atmospheric distances. Over shorter ranges and along paths close 
to the ground, the fluctuations were of the order of 0.1 percent to 10 
percent. 
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The frequency of these fluctuations should be roughly propor- 
tioned to the speed with which atmospheric inhomogeneities are pulled 
broadside through the laser beam or through the field of view of the 
detecting system. The range of energy-bearing frequencies for the laser 
experiments extends from approximately 2 cps to 300 cps for a wind speed 
of 2 to 4 m/sec (figure 15). For given atmospheric conditions, the upper 
frequency limit will depend on the beam diameter or on the field of view 
of the detecting system in a general case. If we take on orbital scan- 
ning speed of 8000 m/sec and a field of view equivalent to 50 meters, 
then frequencies up to 160 cycles per second are to be expected. In 
the computations of detector noise, a figure of 1000 cps was taken as 
the upper frequency limit for the orbital case. 

Other evidence of fluctuation levels to be expected in atmos- 
pheric crossed-beara experiments is available from measured temperature 
fluctuations. Since in the atmosphere velocities vary much less than 
the speed of sound, density variations arise from temperature variations. 
Typically, variations of the order of a few tenths to over one degree 
centigrade are observed; therefore, these variations correspond to 
density variations of up to 1 percent. The number density of aerosol 
particles might be expected to show similar fluctuations. 

In fact, as noted above, there is practically no information 
available on the fluctuations in scattered sunlight recorded by a detec- 
tor on the ground or in a satellite; it is for this reason that we are 
now conducting such experiments. The information obtained will serve 
as input to determine the optimum design for the experiments described 
in later sections of this paper. 

The modulations of radiated power which are produced at higher 
altitudes are equally difficult to estimate. For scattered sunlight, 
variations in aerosol concentrations are likely to be the main source of 
fluctuations at low altitudes. At higher altitudes, fluctuations in air 
density are more important. This is because above 5 km, the aerosol 
concentration decreases much more rapidly with altitude than does the 
air density (figure 16). The contribution to the fluctuating signal, 
for scattered sunlight, is therefore expected to decrease quite rapidly 
with altitude; this will limit the maximum altitude for which the cross- 
correlation technique will permit the signal to be extracted from the 
noise in reasonable integration times. 

This difficulty is not present in the case of ozone emission 
because of the distribution of ozone in the atmosphere. The altitude 
of maximum ozone concentration varies from about 10 to 30 kilometers 
depending on latitude, time of year, etc. Thus, ozone is likely to be 
useful in tracing atmospheric motions at higher altitudes. For the cal- 
culations in this paper, fluctuations in ozone concentration are assumed 
to be similar to those in air density, that is, of the order of 1 percent. 
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In our estimate of the expected modulation amplitude ij^g, we tried to 
account for the decrease of tracer concentration with altitude by choos- 
ing iRMs/I = 10" 3 as a reference value. The associated meteorological 
modulations of the mean radiative power are summarized in table 1. 

3. Detector Noise 


Although the nature of the cross-correlation technique permits 
signals that are buried in the noise to be measured, it is very desirable 
that the above meteorological modulations of the incident radiative power 
produce output signal fluctuations, ij^g, larger than the detector noise, 
i^. In this case, the feasibility of the crossed-beam technique, as 
applied to the study of atmospheric problems, will be determined by the 
magnitude of the background fluctuations and the turbulence length scales 
along each of the crossed beams. The space-time variations of scattered 
sunlight and thermal radiation represent an output from integrated light 
sources and should not vary as much as the local emission and extinction 
processes of interest. The number of correlation volumes along the 
optical path should be small enough to avoid cancellation due to inte- 
grating over too many statistically independent parts. Since turbulence 
is mostly restricted to a few layers, cancellation should be no problem. 
Even if cancellation occurs, it can be partially eliminated by acting on 
the optical path length. Therefore, the main limiting factor appears to 
be the detector noise. 

The detector noise has been calculated for detectors operating 
in the visible (A = 0.55jj,) and infrared (A = 9.6^); the results are sum- 
marized in table 4. 

The sizes of the collecting mirrors are chosen with regard to 
the final signal/noise ratios and are no larger than the mirror systems 
in the 0A0 series of orbiting astronomical observatories. The choice 
of field of view of the optical systems is based on assumed average 
integral length scale of 300 meters. This latter figure agrees very 
well with the scales of clear-air turbulence measured by Burns and 
Rider [21]. To resolve 300-meter scales, a field of view corresponding 
to approximately 50 meters is desirable; this figure at the tropopause 
gives the field of view figures given in the table. 

The values of mean radiance used in the table and discussed in 
previous sections of this paper, together with the field of view and the 
size of the collecting mirror, permit the radiation received by the 
detector to be calculated. 

In the visible region of the spectrum, phototubes or photo- 
multipliers are the appropriate detectors to use, and for the light 
levels in question, the main source of noise in the output signal will 
result from the randomness of the emission of electrons from the photo- 
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Table 4 

Detector Operation and Noise 


Operation 

Visible Light 
(A = 0.55(a) 

Infrared 
(A = 9 . 6 pi) 

Power 

Ground 

500km Orbit 

Ground 

500km Orbit 

Mirror Cross- 
section [cm 2 ] 

10 2 

10 4 

10 4 

10 4 

Field of View 
Minutes of arc (£a) 

4 jt sin 2 (£a/4) ster. 

15' 

1.5 • 10' 5 

0.34' 

7.8 • 10 -9 

15' 

1.3 • lO -5 

0.34' 

7.8 • 10" 9 

Mean Radiance 
[watts/ cm 2 s ter ] 

4.5 • 10- 3 

4.2 • 10" 3 

7.8 • 10 -5 

4 • 10 -3 

Minimum of Mean 
Radiative 
Power I [watts] 

5.8 • 10- 6 

3.3 • 10 -7 

io - 5 

3.1 • 10" 7 

Photon Current 

•—- [Photons/sec] 
hv 

16 • 10 12 

8.9 • 10 11 

4.8 • 10 14 

1.5 • 10 13 

Modulation Band- 
width B [cps] 

300 

1000 

300 

1000 

1 

Detector 

i 

Mean Signal 

L 

Output 
2.6 • 10 -7 
Ampere 

Output 
1.4 • 10" 8 
Ampere 

Input 

10“ 5 

Watt 

Input 

3.1 • 10 -7 1 

Watt 

Detector Noise 

Output 
5 * 10~ 12 
Ampere 

Output 
2.1 • 10 -12 
Ampere 

Noise equiva- 
lent power 
1.7 • 10 -1 ° 
Watt 

Noise equiva- 
lent power 
3.2 • 10“ 10 
Watt 

Atmospheric 

Modulation 

i RMS 

Output 
2.6 • 10 -1 ° 
Ampere 

Output 

1.4 • 10 -11 
Ampere 

Input 

10 _s 

Watt 

Input 

3.1 • 10 -1 ° 
Watt 

Signal/Noise 

’-RMS^d 

52 

6 

59 

1 
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cathode, or "shot-noise. 11 The number of photons incident on the photo- 
cathode may be simply calculated by dividing the energy per photon into 
the incident radiant flux. At 0.55p, the energy per photon is given by 

hv = 3.9 x 1(T 19 watt-seconds. 

If photocathodes of the tri-alkali type are used, the average quantum 
efficiency Q over the visible spectral region will be of the order of 
10 percent. That is, for each ten incident photons, one photoelectron 
will, on the average, be emitted. The cathode current is therefore 
given by 


I d = Q N e e, (29) 

where Q is the quantum efficiency, N e is the number of photons incident 
per second, and e is the electronic charge. The mean square noise 
associated with this cathode current is 


= 2e IB = 32 x lO" 20 I d B, (30) 

where B (cps) is the electronic bandwidth of the detecting system. 

For infrared detectors, the noise is listed in terms of the 
equivalent radiative input power i d , which would produce the same RMS 
value of the detector output in an imaginary experiment, conducted with 
a "perfect" noise-free detector. The best detector to use in the wave 
length range 8 - 12(j, is a mercury-doped germanium- type of photoconduc- 
tive detector [22], which for a detector 0.1 cm^ will have a noise 
equivalent input equal to 


i^(watt) = 1 x lO” 11 n/b cps 1 


(31) 


The botton line of table 4 gives the resulting s ignal- to-noise 
ratio (iRMs/^d) for t * ie f° ur systems, where the modulation in the 
incident intensity at the detector due to the atmosphere is assumed to 
be 0.1 percent in each case. It appears that the feasibility of the 
application of the crossed-beam technique to the atmosphere will not 
be limited by detector noise unless the signal fluctuations caused by 
atmospheric effects are much smaller than presupposed here. 

VI. THEORY OF CROSSED-BEAM SCANNING 

In the global coverage of vast and sparsely instrumented areas, 
environmental monitoring systems would have to gather information 
quickly enough that the results can be used for the prediction of 
weather and the prediction of navigational and communication hazards 
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such as clear-air turbulence. This consideration leads to the exten- 
sion of the crossed-beam concept from space-fixed beams to scanning or 
sweeping beams. M. Wolf [2] has already successfully operated a scan- 
ning crossed-beam system. He rotated two ground detectors in such a 
way that the beam intersection point moved along a horizontal path at 
approximately 100 km altitude. A cross correlation of the optical 
fluctuations from the night-glow emission gave the height of the air- 
glow phenomenon with great accuracy. We now extend the idea of crossed- 
beam scanning to include systematic beam separations in space and time 
and the measurement of turbulence parameters. We will indicate how 
the two-beam product mean values from moving beams are to be inter- 
preted in terms of convection speeds, turbulence length scales, and 
space correlations (wave number components). The theoretical treatment 
is kept general and could be applied to orbital as well as ground 
detectors . 

1 . Geometry of Scanning Beams 

We define the scanning or sweeping motion of crossed beams by 
the motion of the points of minimum beam separation A and B. The 
straight line which goes through these two points had been introduced 
as the beam normal. In case of crossed-beam scanning, this normal 
moves broadside and sweeps out a certain surface. We have chosen to 
classify crossed-beam systems in terms of the motion of the beam normal 
and through the orientation of the beams relative to the surface of the 
sweeping normal. This gives a unified treatment which can be applied 
to both orbital and ground detectors. The following discussions are 
restricted to a crossed-beam system where both the speed of the normal 
and the orientation of the beams relative to the surface of the normal 
are time invariant. Such crossed-beam systems, shown in figure 17, may 
be compared to the translation of a rigid mechanical structure. The 
beam arrangement is thus uniquely classified in terms of the following 
invariant parameters: 

(a) The scanning speed U-g describing the speed of the 
point B of minimum beam separation inside the sur- 
face of the sweeping normal. 

(b) The inclination (3 of the path of B relative to the 
beam front, that is, relative to a plane which is 
parallel to both beams. The path of point B will 
be called the "intersection path." 

(c) The inclinations and a-g of the two beams relative 
to the surface of the sweeping normal. 
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The motion of the individual beam is now described in two dif- 
ferent space-fixed coordinate systems, also shown in figure 17. The 
first system is needed to describe the position of a beam point rela- 
tive to the earth’s surface. The z-axis points vertically upward and 
the x-axis follows the line of minimum beam separation. The y-axis 
follows from the condition that the x,y,z system should be an ortho- 
gonal system. The coordinates of the second system are defined 

in terms of the position which the sweeping beams occupied at the 
instant t = 0, that is, at the start of the experiment. The |-axis 
follows the beam normal, the rj-axis the leading beam, and the £-axis 
the trailing beam. This system is not necessarily orthogonal, since 
the angle 


a AB = ' °A ' °B (J2) 

between the beams is not necessarily a right angle. The relation 
between the two coordinate systems follows from geometrical considera- 
tions and may be given as follows: 

^(f) = + I 

y(I) = y Q + ti cos - £ cos a A (33) 

z(T) = Z c + T) sin Og + £ sin ct^. 

Here x Q denotes the origin of the system, that is, the position 

of the moving point B at time zero. The trajectory of point B follows 
the "intersection path." Each point on this path shall be denoted by 
the path length p relative to the point 5? 0 . 


*B (P > 



x o 


sin p 

= 

y 0 

+ p 

cos (3 


L z °_ 


0 


(34) 


The scanning speed Ug describes a velocity vector which follows the 
intersection path. This vector may thus be denoted by 



u s 


sin p 

= 

V s 

= U B 

cos p 


W s 


0 


(35) 
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Similar trajectories can be written for any other point of the trailing 
beam "l". At time zero this beam coincides with the £-axis and may 
thus be found at position x(0,0,£). The trajectory of this point is 
then given by 





x D + p sin p 

X^.p) = x(0,0,g) + ^[p(t>] = 

yi(^t) 

= 

y Q - £ cos a A + P cos p 




z Q + g sin a A 


(36) 

The fluctuations i-^t), which are recorded for the trailing beam, follow 
from the optical integration along this beam.* 


ii(t) -<I X > / k(x!(^,t); t) d£. 


(37) 


Similarly, an arbitrary point on the leading beam ff 2 M is denoted by 
the beam length q and moves along the trajectory. 


X2(|,n,t) = x(g,T),0) + ? B [p(t)] = 


x 2 (5,T],P) 


x Q + | + P sin p 

yaCM.P) 

= 

y Q + T] cos a - £ cos « A + P cos p 

Z2(i, Tj,p) 


z D + r] sin Og + t, sin 


(38) 


The optical integration along the leading beam can thus be expressed by 


=<i 2 > / k ( x s(M, p); fc ) 


(39) 


Light source fluctuations are neglected since they are mostly uncorre- 
lated; i.e., a further noise contribution, which is suppressed by cor- 
relation techniques. The operator < > is used for the time average of 
moving beams. 
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The two-beam product mean value of a sweeping crossed-beam system follows 
from the time average of the instantaneous product 

00 

R(i> t = 0) = <ii(t) i 2 (t) = J jT<k[x 1 (^,p(t)) ; t] k [x 2 (|, t], p(t) ) ; t]>d^dr] 

-00 

00 

= //vB ( 4 °) 


One finds that the two-beam product mean value is an area integral over 
a great number of observer pairs each of which travels parallel to the 
intersection path. 

2. Length Averages of Frozen Patterns 

The individual observer averages over the fluctuations which he 
encountered along his trajectory. The fluctuations are produced by scan- 
ning inhomogeneous spatial distribution of radiative power and by true 
temporal variations as communicated by a space-fixed observer. We have 
thus a mixed space and time average. This averaging process is analogous 
to the way the human eye evaluates a TV picture. Each frame of the 
television is displayed on the screen by a rapidly moving electron beam. 
To the viewer the modulations of the electron beam are displayed so 
rapidly that the resulting image appears as a "frozen pattern"; that is, 
like a flashlight image. Although the information which is displayed by 
the scanning beam was recorded as a time function, its interpretation 
assumes a space history. Only by comparing several frames is it possible 
to detect generation, convection, and decay of new patterns. 

We now assume that an observer pair scans so rapidly that the 
associated two-point product mean value closely approximates a length 
average R? along the observer trajectory. This means replacing the time 
history of a moving beam point with the flashlight image (t = 0) at all 
points P along its trajectory. 

=<k[x 1 (^,p(t)>; t] k[x 2 (|,Ti,p(t)); t] > 



P oo 


p=P 

r 

j 

p=0 


Mx^.p); t = 0) k(x 2 (|,r),p); t = 0) dp. 


(41) 


For convenience, this length average will be denoted by the operator 
"-p", i.e., ~ Rj?. The last approximation is valid whenever the spatial 
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distribution or "pattern" of k is frozen [k(i?, t) = k^O)] over the time 
period of the scan. It is also valid if a frozen pattern is convected 
at a constant speed Uc, since this would alter only the effective scan- 
ning speed of the moving observer pair. The last equation is therefore 
valid for all convected frozen patterns; that is, for all possible direc- 
tions of the convection speed. Furthermore, it can be shown that the 
pattern does not have to be entirely frozen. The approximation of length 
averages is still good as long as the temporal variations which are 
experienced by a third observer travelling with the convection speed can 
be treated as small perturbations and are statistically stationary. 

Replacing the time average of scanning observers with a length 
average over the observer's trajectory leads to the following expression 
for the two-beam product mean value of a rapidly scanning system: 

00 

<ii(t) i 2 (t)> s =<I 1 XI 2 > '/ R£(?) drj d £ = 

-00 

(42) 

00 

<Ii><I 2 > 77 t = 0) k(x 2 (i,q,p); t = 0)) dq d£. 


In the above averaging procedure, it was always assumed that no 
time lag was introduced in the data reduction. However, in space-fixed 
crossed-beam arrangement, the systematic variation of time lags was 
required to separate convection and generation and decay processes and 
to measure the convection speed. Similar interpretations of temporal 
correlation, which hold for moving beams, will now be discussed. 

The fluctuations of radiative power which produce the signal 
i 2 (t) along the leading beam have been on another position imag at 
the time t - t. In case of a convected frozen pattern, this earlier 
position is parallel to the leading beam and may be found by tracing 
back the scanning motion Us and the convective motion tfj-, — (Uq, Vq, Wq) . 
The signal of the delayed beam i 2 (t “ t) is thus identical with the 
signal from an imaginary undelayed beam, each point of which could be 
found by the translation 


*2,imag. ( ^ ; P ’ T) = *2<M; t(p)) 


(if + U c ) T. 


(43) 
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It follows that the above optical integration of two-beam product mean 
values, equation (42) can be extended to the analysis of a delayed beam, 
if one replaces x 2 with x 2 , imag* 


<Ii > <Ip> 


k(xi(£;p);°) k(x . ( E; r), p ; t) ; 0) dt]d(;. 


(44) 


The interpretation of this two-beam product mean value follows in analogy 
to the discussion of space-fixed beam by considering the two-point product 
mean value between the two points Aj_ and which are defined through 
the minimum separation |(t) between the trailing beam and the imaginary 
beam. The position of these two points is indicated on figure 18 by the 
projections of the imaginary beam into the side view and top view of the 
plane of the sweeping normal. The figure shows that the two pairs of 
points A^, Bj_ and A, B differ by a translation and a change of separation 
distance. The translation occurs inside the beam front and may be charac- 
terized by its Z and Y components. 


(V s + V c ) T 

Az(-t) - z Bi Z B cot a A + cot 


(45) 


(V s + V COt °B 

Ay(-r) = y B . - y B = cot ^ + cot ' ^ 


(46) 


The change in separation distance is given by 


U S +U C 

AS(-x) ^Ai " V + v r ^ 

S t 


“b 


cot a. + cot 
A 


— (U s +U C )(-T). (47) 


As long as the frozen pattern of k is statistically homogeneous over the 
beam front, a translation by -Az, -Ay would not affect the two-point 
product mean value. However, this translation brings the pair A^B.^ back 
to the original normal (i.e., r\ = £ = 0). In fact, the translated point 
Bi coincides with B. The translated point A^, which is denoted in fig- 
ure 18 by C, defines the new beam separation 


£*(-t) = £ + A| (- t) = £ + C AB (U B sin p + U c )(-t), 


(48) 
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where the constant abbreviates the cotangent factor 


cot 

C AB cot a A + cot Og 


( 49 ) 


The only effect of a time delay is ? thus to change the actual beam 
separation, g, to the effective, The approximation of two-point 

product mean values may therefore be extended from space-fixed to moving 
beams in the following way: 


<i 1 (t) i a (t + t)>^ r£ T| = 0; £ = 0) 

<i 1 (t) i 2 (t)> | ~ r£ a; I] = °, £ = 0) 


(50) 

k(x Q + p sin p, y 0 +P cos p; z Q ; 0) k(x 0 + p sin P + |*(t);Y 0 +P cos p, z Q ; 0) P 
k(x Q +p sin p, y 0 + p cos p; z o ,0) k(x Q +p sin p + £; x Q + p cos p; z Q ;~ 0) P 


In a physical sense the time lag x M is needed for the combined 
motions U s + Uc to cover the original beam separation, g . The delayed 
leading beam then intersects the undelayed trailing beam and the corre- 
lation becomes a maximum, since the two points of minimum separation 
sweeps the same eddies. The convection speed is thus measured by 
averaging over a row of eddies along the intersection path whereas 
space-fixed beams average over a single row of moving eddies, which 
pass a fixed observation point at different times. 

3 . Interpretation of Two-Beam Product Mean Values 

The above approximation of space correlation curves, equation 
(50), allows several important applications. The first application is 
the measurement of convection speeds. It follows from the fact that 
space correlation curves have a maximum at zero space lag. Therefore, 
the time lag T^(g), which is indicated by the maximum of the temporal 
beam correlation curve<i 1 (t) i 2 (t + t)^ 5 should correspond to g* = 0. 

E* V ■ 0 - 5 + C AB ( °C + U B sin » V <5l) 
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This is an equation which may be solved for the group velocity component: 


“ c ' ‘ “b sin p - <52) 

The associated estimates of group velocities are graphically illustrated 
in figure 19. In equation (52), the second term, Ug sin p, corrects for 
the beam motion. If this term is large compared to the group velocity, 


I Ug sin p/U c | = |U g /U c | g 5, 


(53) 


then the measurement of Uq becomes difficult since it involves the sub- 
traction of two large quantities of about equal magnitude. The scanning 
speed Ug must be much larger than the convection speed to justify the 
assumption of a frozen pattern. Therefore, the last inequality implies 
that p must always be a small angle. In other words, the beam front 
has to move almost parallel to itself. 

The second application is the estimate of the integral length 
scale L x . Since the two-beam product mean value is proportional to the 
space correlations curve RP(£"), the normalized space correlation curve 
can be approximated by normalizing the two-beam product mean values with 
their maximum 




1 V£\ 

u c ) 


Rg(£*,0,0) 

(o,o,o) 


<ii(t) i 2 (t + t) >£ 

9 


<ii(t) Mt + 


(54) 


and by plotting this coefficient as a function of the effective space 
lag Once the convection speed Uq is known, this lag may be calcu- 

lated without using the scanning speed, Ug. Substituting equation (51) 
into equation (48) shows 


i* = 


T M ( ^ 

Un 


(55) 
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The plot of space correlation functions is illustrated in figure 20. 
The integral length scale, L x , in the direction of the beam normal 
follows from integration 


L = 77 


0(0 df 


( 56 ) 


The third application is the estimate of three-dimensional wave 
number components. The components result, if one wants to describe the 
space history of a frozen pattern of extinction coefficients k(x > , t = 0) 
through a superposition of harmonic waves. Each of these waves is 
characterized by the number of cycles per unit length \|/ x , \Jr 2 , which 
may be oriented in either the |, rj, or £ direction. The individual 
wave number component, S^., denotes the mean square amplitude of such an 
individual wave and follows from a three-dimensional Fourier transform 
of the two point product mean value. 


00 

S k ( ^ Iff d£dr)di. 

— 00 


(57) 


For wave fronts which are parallel to the plane of the beams, we get 
^2 = ^3 = and the inner two integrals describe nothing but an inte- 
gration over the beam front which is automatically performed. 




^ dr)d£ 



— 00 


R £(?* ) dr]d£ 


<ii(t) i 2 (t + t)>| 
<I ± > < I 2 > 


e 


-mi. 


(58) 


Substituting into the previous equation gives the three-dimensional wave 
number component by a one-dimensional integration 


V+i.o.o) ■ 5js 


<il(t) i 2 (t + 

I-i Ip 




■>v 


d | 


ic 


(59) 
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All properties of the integrand are experimentally accessible 
to a moving crossed-beam system. 


The fourth application is an estimate of the normalized power 
spectrum. Consider several crossed-beam systems with different beam 
separations, | | l5 | 2 , ... , which scan the same integration path 

simultaneously. According to equation (54), the associated temporal 
correlation curves R s (x)| = < i x (t) i 2 (t + t) >|^ should all collapse 

into a single curve if they are translated along the time-lag axis until 
the maximum coincides with the ordinate. In other words, each of the 

individual space- time correlations R s (t)£ could be generated by trans- 

5 i 


lating the autocorrelation R s ( T )| = o t * le amount t^(|). This is 

indicated in figure 19 by the dashed curves and follows from the assump- 
tions of a frozen pattern. Therefore, by replotting the normalized two- 
beam product mean value of a single run | as a function of the 

"retarded time lag," 


T* = T - Ijjdp (60) 

should give a curve 0(t“) which closely resembles the normalized auto- 
correlation function. The Fourier transformation of this curve should 
then give an estimate of the normalized power spectrum 


S k ( f ) 



-i2jtfr* , * 
e dr . 


(61) 


The consideration of several simultaneously estimated time correlations 
curves R s ( T Hi also implies a lower limit on the scanning speed. The 


assumption of a frozen pattern is obviously wrong if the individual 
curves R s (x)| depart considerably from the translated autocorrelation 

R s( T )|=o* This case is illustrated in figure 19 for | = I 2 • The con- 
cept of a frozen pattern therefore seems justified only as long as the 
range of time lags used is small compared to the eddy lifetime. In the 
atmosphere, eddy lifetimes are usually in excess of 5 minutes, except 
for layers very close to the ground. On the other hand, our wind tunnel 
tests suggest that one would probably have to scan over 50 statistically 
independent regions to retrieve the local information. An average 
integral length scale of Ly=300m [21] gives a region of L = 50L y = 15 km 
length. To scan this region in a tenth of the average lifetime, one 
needs a scanning speed in excess of 15 km/0.5 min = 0.5 km/sec. 
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4. Measurement of Altitude Profiles 


It remains to consider the altitude variation that is introduced 
when. setting a time lag. The convection speed is determined from a 
single time lag and its value belongs to the altitude that is indi- 
cated by t m , 


V s + V C 

h(T M ) = h (T = 0) + co't coT^' V 


(45) 


If one wants to measure an altitude profile of convection speeds 
in a single sweep, one needs several independent detectors which view the 
altitude range simultaneously. Each of these detectors should indicate 
a correlation maximum at a different One easy solution is suggested 

by equation (51). All it takes is to set the photodetectors at different 
inclinations (3 against the flight path. Since all of these inclinations 
are small, as discussed previously, this can be accomplished using only 
one telescope. One simply collects the radiative power from several 
off-axis pin-holes in the focal plane of the collecting mirror, as shown 
in figure 21. 

Multiple recorders are mandatory when measuring the space cor- 
relation R s (t*). The calculation of integral length scales, equation 
(51), and wave number components, equation (59), requires a knowledge 
of the spatial beam correlation R s (|*) over the entire range of space 
lags, 



(62) 


where 0(| w ) is not zero. The associated range of time lags is fixed by 
the convection speed and cannot be shortened by special test arrangements. 


T 

max 


T . 
mm 




X 


C 


in 



(63) 


Typical time lag intervals may be derived setting - Imin < - 1000M 

and Uq ^ lOOM/sec. They are of the order of 10 seconds, which is much too 
long when considering the altitudes that are covered meanwhile by scanning 
speeds in excess of 0.5 km/sec (orbital speed 8 km/sec). Instead, one 
has to break up the required time lag interval into subintervals of 
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length At, which are derived from the permissible altitude error £h as 
given by equation (45). 


cot a. + cot QL 

At = — — (64) 

s C 

One independent photodetector $, = 1, 2, 3, is required to cover the 

lag interval - 1) At S t § f A t. The inclination (3g of the individual 
detector has to be chosen such that a maximum of the time correlation 
R s (x) is falling into the interval. This leads to the condition 


T M (Pg) » U - \) At 


(65) 


which is sufficient to calculate Substituting equation (65) into 

equation (51), solving for (3^ and using the expressions for equa- 

tion (49), and At, equation (64), gives the expression 


sin 



+ (v s + v c ) tg % 



( 66 ) 


Choosing the inclinations of the individual photodetectors according to 
equation (66) guarantees that a maximum t h(P^) will be found (that is, 
a group velocity measurement) in atmospheric layers of thickness At. 
Furthermore, the associated retarded space lags £* = £C*M((3g)) are 
spread throughout the desired range of retarded space lags. This pro- 
cedure also reduces the requirement of homogeneity inside the wave front. 
Homogeneous behavior is demanded only on the now small translations 
(At) and Ay (At) . 


VII. STEP-BY-STEP SIMULATION OF IN-FLIGHT EXPERIMENTS 


Having developed the interpretation of scanning beams, it remains 
to indicate how the underlying motion of a rigid beam structure can be 
activated in practice. The measurement of convection speeds, turbulence 
length scales, and wave number components requires a beam normal which 
moves broadside through the atmosphere at a scanning speed in excess of 
0.5 km/sec. Furthermore, the telescope which defines this normal was 
restricted to a translation. However, a rotation of the beam structure 
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around the normal will not change the optical integration over the beam 
front. Therefore, the telescope motion may be characterized by the con- 
ditions C^b, and P must be invariant. Such motion can be achieved with 
both orbiting and ground detectors. The associated test arrangements 
are illustrated in figure 22. Using such arrangements has the great 
advantage that orbital experiments can be simulated and checked out in 
ground experiments . 

1 . Test Arrangements 

For orbiting telescopes, the beam separation distance, g, can 
be kept constant if the two detectors orbit at the same speed, that is, 
at the same height. The path of the normal is then horizontal and the 
surface of reference coincides with surfaces of constant altitude. For 
ground detectors, a very similar sweep can be arranged by separating the 
two telescopes by a baseline b and by rotating both beams at the same 
rate (QA = Og = a = const). In this case, the normal is also horizontal 
and sweeps along a circular arc such that the surface of reference almost 
coincides with a horizontal plane. The main difference between orbiting 
and ground detectors is not in the individual sweep but in the way in 
which the length averages of the two-beam product mean value are estab- 
lished. Orbiting detectors allow averaging over a single horizontal path 
of infinite length. Ground detectors average over a finite portion of a 
circular arc, the length of which is determined primarily by the telescope 
baseline b. However, infinite statistical information can still be col- 
lected if one sweeps repeatedly and includes these repeated runs in the 
averaging procedure. Both averaging procedures give identical results 
if the pattern of turbulence is spatially homogeneous over the interval 
of the orbital sweep or if it can be made homogeneous by subtracting 
large scale, inhomogeneous "means." 

2. Preliminary Specifications 

We have tried to specify the crossed-beam test arrangements, 
shown in figure 22, in terms of the independent variables which describe 
the telescope operation and/or adjustment. These variables are as 
follows : 


A optical wavelength set by monochromator 

b baseline or distance between telescopes 

O' beam intersection angle 

AB 

scanning speed that is the speed of the moving 
normal 

| beam separation along beam normal defined by points 

of minimum beam separation 
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T electronically introduced time lag between beams 

(3 inclination between the plane of the beams and the 
path, which is followed by the point of minimum 
separation 

N number of desired altitude positions, defining number 
of photodetectors in local plane of telescope 

P interval of intersection path used in averaging 
procedures 


Table 5 


Independent 

Variable 

Minimum or 
Lower Limit 

Maximum or 
Upper Limit ; 

Tolerance 
Given by 

Comment 

Optical 
Wave length 

0.38n 

0. 72p 

Entire window 
for particulate 
tracers 

Visible window 

A 


12p 

For gaseous 
tracers 

Infrared window 

Telescope 

Separation 

b 

2 cot a h . 

min 

2 cot a h 

max 

Position error 

A 0 = Ac 

h is the distance 
between the beam 
intersection path 
and the baseline. 

Beam Intersection 
Angle o- AB 

15° 

75° 

Altitude error 

o Ah 

/jo: - 2 cos'-a^ 

b 

Depends on point- 
ing accuracy 

Beam Separation 

- 2 Lx 

21, x 

A| c % 1^/ 5 
(beam diameter) 

is the length 
scale. The beam 
diameter factor 
1/5 is obtained 
from wind tunnel 
experiments . 

Time lag 

0 

^max ^mix/ U C 

“ 2L x/U c 

At = Lx/5U B 

L x and Ug depend 
on applications . 

Inclination of 
intersection 
path with beam 
plane |p| 


10 k 

Ag - 2A| 0 /b 
= 1 L x 

Equation for 

C AB = sets 

angle of view or 

horizontal point- 
ing accuracy. 

Number of 
photodetectors 

5 

depends on 
hardware and 
receiver 

- 

At least five 
points on a curve 
to find a maximum 

Length of 
intersection 
path p 

Not known 
length scale 
of homogeneous 
trend > 1^ 

Not known ~ 
depends on 
correlation 
maximum and 
receiver noise 

^ * L x 

One sweep from 
orbit or pieces 
of repeated sweeps 
from ground 

Scanning 

speed 

0.5 km/sec 

8 km/sec 
(Orbital speed) 


Tolerance depends 
on application 
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The desired test arrangement can now be specified in terms of 
the lower limit (minimum), upper limit (maximum), and tolerance (A) 
which are required for each of the independent variables. These limits 
and tolerances depend on the desired applications. The considerations 
of the previous sections allow us to establish rough estimates, the 
expressions of which are listed in table 5 . 

The minimum spacing of space lags, A| 0 , has been determined in 
our wind tunnel tests, which seem to indicate that A| Q = 1/5 L x gives 
satisfactory space correlation curves. The tolerance on time lag, At, 
is identical with the sampling period of the A/D converter. It has been 
approximated by the time the beams would need to scan over a space lag 
interval . 


L 

Ax = A| 0 /U b « . (67) 

B 

Some specifications for the collector optics (telescope) follow 
from the tolerance column. The angle of view or pointing accuracy is 
given by the smallest value of or Aft . The rotation rates required 
for stabilization of the telescope platform follow from the condition 
that the beam oscillation should be kept so small that it takes more 
than the sampling period. At, to produce an angular error of size 

jW or To" 


10 AT 


2 A U B 
cos^a -r- ~ 

b L x 


Ap 


g -L Ag 

10 At 


jl!b 

10 b 


( 68 ) 

(69) 


Typical pointing and tracking accuracies have been calculated 
for orbital and ground-based detection (see table 6). 

The above considerations are sufficient to outline telescope 
arrangements needed to scan a certain altitude from either the ground 
or orbit. However, the characteristics of the crossed-beam system in 
terms of the length, P, of the scanned portion, the number of photo- 
detectors, N, and the spatial and temporal resolutions of convection 
speeds, length scales, and wave number components cannot be determined 
without conducting several experiments. 
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Table 6 


Pointing Accuracies and Stabil ication of Telescope Platform 


Pointing 


AP ^ | -p KXfk 2 cos 2 0 ! y' 


S tab i 1 ization 


2U. 




B . 
10b 5 


z^a ^ cos^a 


^ U R 

2,- B 


b U, 


Detector 
Pos it ion 

Base- 
line 
[kms ] 

Turbulence 

Length 

Scale 

L x [m] 

Altitude 

Tolerance 

Ah[m] 

Accuracy 

ground 

10 

100 

50 

= Af3 = 4 x 10 -3 = 
15 arc min 

(U B = 500 
m/sec) 




AD! = Ap = 10“ 2 sec -1 = 
34 arc min/sec 

500 km 
orbit 

1000 

nrn 

i on 

Aa = A£ = 10" 4 = 0.34 

arc min 

(U B = 8000 
m/sec) 

250 

1UU 

Aa = AP = 1.6 • 10 -4 sec -1 
0.5 arc min / sec 


Comment: 


a a A “b 2 a AB 


45°. 
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We plan to conduct several experiments in which the true range 
and resolution of scanning crossed-beam systems are established by com- 
paring the optically measured convection speeds with horizontal wind 
measurements from meteorological towers and radiosondes. These experi- 
ments represent a step-by-step approach from simple space-fixed beams to 
rapidly scanning beams. 

3. First Generation Experiments 

The two sensors of the crossed-beam arrangement are placed on 
a short baseline with their intersection point set at a moderate height 
(about 25 m) above ground so that the results of the experiments can be 
compared to tower-mounted anemometers. The intersection points are fixed 
in height, but the beam separation distance is varied from 0 to 200 meters 
in a number of discrete steps. At each step the following results are 
obtained : 


(1) The rms intensity of the registered fluctuations. 

(2) The detector s ignal- to-noise ratio. 

(3) The horizontal wind speed, turbulent lengths, and 
wave length spectrum of the phenomena occurring in 
the region of beam crossing. 

(4) The spatial and temporal resolution of the system. 

(5) The effect, if any, of the beam diameter on the 
resolution of the system. 

The first series of tests should serve to establish the sound- 
ness of the apparatus design. Once satisfactory results are obtained 
with this short range arrangement, the baseline of the detectors will 
be increased to approximately 20 km to produce intersection heights of 
up to 10 km. The beam intersection point will be adjusted to discrete 
altitudes above tower heights to assess (1) rms levels and frequency 
range of fluctuations, which originate at higher altitudes, (2) the 
influence of the variation of the beam intersection angle on the meas- 
ured correlation values, and (3) the scale of statistically inhomogeneous 
fluctuations . 

These experiments will sense phenomena over longer ranges and 
under conditions of reduced light intensity and tracer concentration and 
will therefore be a more severe test of the feasibility of the technique 
as applied to atmospheric phenomena. The beam separation distance, g, 
should be varied to approximately 600 meters (approximately twice the 
mean eddy sizes at an altitude of 10 km) and data will be obtained on 
the five parameters listed at the beginning of section VII. 2. 
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The considerations of section V anticipate no signal detection 
problems up to about 15 km. If the measurements are to be extended to 
an altitude of 30 km, a point might be reached where either the tracers 
or the light intensity or both fall below the minimum value necessary 
for producing a detectable signal. Since the ozone concentration at 
a given height varies with weather conditions, only actual field experi- 
ments will define the height limitations of the system. 

4. Second Generation Experiments 

Once the feasibility of applying the crossed-beam technique to 
the measurement of atmospheric turbulence phenomena has been established 
by the series of tests using fixed beams, an attempt will be made to 
obtain continuous coverage of turbulence phenomena by scanning the beams 
over the region of interest. 

The two beams will be scanned in elevation so that their inter- 
section point will describe horizontal traces in the region of interest. 
The scanning should require less than 1.5 minutes so that the turbulence 
phenomena can be considered as frozen during the period of one scan. 

The main point of interest in this series of tests will be to 
determine the effect of the scanning speed of the beams on the measured 
quantities. In particular, the replacement of time averages in a fixed 
spot with length average over the intersection path has to be tested. 
Certain types of inhomogeneous atmospheric turbulence are conceivable, 
where the average over one long path interval may have to be replaced 
with an ensemble average over many repeated sweeps of short length. 

The occurrence of such turbulent fields and the difference between a 
single sweep and multiple sweeps can be explored with rapidly scanning 
ground detectors. Multiple sweeps with orbiting detectors will not be 
possible. 


A single sweep with two orbiting detectors will be sufficient, 
as long as the inhomogeneity is restricted to large scale phenomena. 

The inhomogeneity can then be removed as a nonstationary mean [23] or 
by employing structure functions [24]. Available experimental evidence 
[24] suggests that the suppression of nonstationary trends will be pos- 
sible for most atmospheric applications of practical interest. 

A second point of utmost practical importance is a decision 
whether an orbital sweep with two satellites could be replaced with 
the flight of a single satellite carrying two telescopes, A and B. 

The possibility of finding a crossing height by triangulation between 
the telescopes of a single satellite is provided by a time lag t in the 
data reduction. This triangulation is illustrated in figure 22. The 
radiation sources which constitute the leading beam B = Ai at time t - t 
could be represented by an undelayed, imaginary satellite which is 
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trailing on the same orbit and possesses one telescope Ai = B. 

The line of sight of this delayed or imaginary beam A^ = B intersects 
the undelayed beam A at an altitude which is given entirely by the time 
lag t and the superposition of scanning and convective motions. Thus, 
altitude profiles can be measured by setting different time lags as 
discussed in section VI, 4. 

The sweep with a single satellite is possible whenever the 
atmospheric pattern remains frozen over the entire range of time lags. 
For = Og = 1/2 qab = 45°, the maximum time lag is identical to the 
time, 2b/UB, the satellite needs to travel twice the altitude. This 
lag is equal to 1000/8 = 125 seconds. Many atmospheric patterns will 
be frozen over this time period. Therefore, a crossed-beam mission with 
two telescopes mounted on the same spacecraft appears feasible. The 
astronomical telescope (ATM), which is built by MSFC, would meet all the 
requirements outlined in this report and could thus be used for crossed- 
beam missions. The difference between astronomical observations and a 
crossed-beam experiment is that two telescopes are needed instead of 
one, and that a new mount is needed to point the telescopes downward 
instead of upward. 


CONCLUSIONS 


The ultimate objective of the flight experiment is to improve 
numerical methods for the prediction of weather, navigation and communi- 
cation hazards, and pollution by providing the following data: (a) 

horizontal winds at global positions and altitudes which are requested 
by weather prediction centers and (b) turbulence spectra and scales 
showing the exchange of kinetic energy between large scale and small 
scale atmospheric motions. Other conceivable objectives, which are not 
treated in this report, would be (a) the early warning of clear-air 
turbulence by sweeping selected flight paths, and (b) water budget sur- 
veys employing a crossed-beam arrangement of millimeter wave detectors 
in this spectral region. All of these objectives may be met on earth 
orbital missions. If these missions should be successful, several 
generations of future missions become conceivable, which would monitor 
the ionosphere, the earth's magnetosphere, planetary atmospheres and the 
outer layers of the sun. 

The crossed-beam method is a new test arrangement of two con- 
ventional remote sensing devices. Triangulation between the lines of 
sight of two radiation detectors is employed to select a local flow 
region which is centered around the line of minimum beam separation. 

Each of the two detectors will react to all emission, scattering, or 
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absorption processes inside its entire narrow field of view and the 
problem is thus to retrieve the local information from the region where 
the two beams cross. To do this the time history of the received radia- 
tive power of each beam is recorded and used as an input for statistical 
data reduction programs, which are commonly used in random vibration 
analysis and communication studies. The resulting time correlation 
between the two signals partially eliminates the integration over the 
optical path and retrieves the local information. The crossed-beam 
method works whenever detectable modifications of radiative power occur 
which are common to both lines of sight. 

The crossed-beam concept has the following potential: (1) The 

art of remote sensing is extended to cover dynamic effects such as 
disturbances, turbulence, and wave-front phenomena in situations where 
other methods are restricted to single-point or single-path experiments 
and average over the optical path. (2) The method compares the fluctua- 
tions at the two points of minimum beam separation in a way required in 
the analytical description of random fields. (3) Crossed-beam spectros- 
copy should give the mean values of thermodynamic properties in inhomo- 
geneous, turbulent, or unstable media which are not accessible to direct 
sensors. (4) A single sweep of frozen atmospheric patterns may estab- 
lish the complete altitude profiles of turbulent intensity, turbulence 
length scales and wave member components. (5) The method is extremely 
versatile, since the tracking of man-made radiation sources is not 
necessary. A passive operation is possible using a natural radiation 
background. Furthermore, the entire electromagnetic spectrum from gamma 
rays to microwaves could be used. 

The crossed-beam method has been experimentally verified in sub- 
sonic jets by comparing optically measured profiles of convection speeds, 
turbulence scales, space-time correlations and spectra with conventional 
hot-wire results [6]. Crossed-beam scanning over atmospheric distances 
has already been successfully demonstrated by measuring the height of 
the airglow phenomenon [2]. The crossed-beam method will work only if 
the following requirements are met: 

(1) The atmospheric components of tracers must produce local 
changes in radiative power which are not lost during transmission to 
detector. In the atmosphere the mean radiative power is supplied by 
extended natural sources such as scattered sunlight and thermal emission. 
Inside the atmospheric window meteorological modulations occur, whenever 
a single beam sweeps over an inhomogeneous distribution of emission or 
extinction centers such as aerosols, ozone or water vapor. 
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(2) Light source fluctuations, shot noise, and instrument noise 
must be reduced to mean square levels which are not several orders of 
magnitude greater than the integrated signal. This is a far less restric- 
tive requirement than the normal, which is that the erroneous noise should 
be at least one, and preferably two or three, orders of magnitude less 
than the genuine signal. The root mean square levels of the above mete- 
orological fluctuations have been estimated and appear to be always 
larger than the noise levels of already available detectors such as the 
astronomical telescope (ATM). 

(3) The lines of sight must intersect at or scan over a 
selected position for a sufficiently long integration period to pull 
the correlated signals out of the light source fluctuations, flow fluc- 
tuation, shot noise, and instrument noise. The necessary path length 
is intimately related to the magnitudes of the various noise contribu- 
tions to the detected signals. Although, in principle, any degree of 
noise can be eliminated with suitable integration periods, the finite 
record lengths and the dynamic interchannel displacement set practical 
limitations to retrieving the local signal from the integrated signal 
and the combined light source noise, shot noise, and instrument noise. 

In a very rough estimate, a flight path interval of 15 km length has to 
be scanned with speeds in excess of 0.5 km/sec to measure horizontal 
wind and length scale profiles. Ground experiments are urgently needed 
to determine signals and noise for various "crossing heights." 

On occasions the problems on meeting the above requirements will be 
formidable and may prohibit some of the desired applications. However, 
since we can pick the detected radiation from any region of the electro- 
magnetic spectrum and since atmospheric applications allow large beam 
diameters and arbitrary orientations of the two lines of sight, we 
believe that there are many experiments where the above problems may be 
overcome, and useful, hitherto unavailable, data obtained. 

The interpretation of the beam correlation (two-beam product mean 
value) is subject to the following theoretical restrictions. 

(1) The number of statistically independent regions and gas 
dynamic phenomena is small enough to avoid cancellation of local fluctua- 
tions which would occur when integrating over too many statistically 
independent parts. 

(2) Either emission or extinction dominates the local change 
of radiative power inside the correlation volume. 
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(3) The integrated fluctuation of the generalized extinction 
coefficient is sufficiently small to permit linearization. However, 
large changes in the mean value of the transmitted radiative power are 
permissible. Also, large local fluctuations are tolerable if one inte- 
grates over several statistically independent regions. 

(4) The interpretation of a sweeping erossed-beam system is 
possible, if the region of interest is scanned rapidly enough to assume 
a frozen pattern and if the emission and/or extinction fluctuations are 
statistically homogeneous over the altitude tolerance or step which is 
desired in the altitude profile. 

(5) The fluctuations of radiative power are locally isotropic 
over a "correlation disc." This assumption is sufficient, but violations 
are often tolerable. Crossed-beara experiments in anisotropic and inhomo- 
geneous jet shear layers here already yielded a good approximation of 
two-point product mean values. 

A step-by-step simulation of an in-flight experiment is outlined, 
which would allow one to assess the effect of the above practical limita- 
tions and theoretical restrictions. This is simulated by rotating two 
ground detectors in such a way that their intersection sweeps along a 
horizontal arc with orbital speed. Preliminary specifications for both 
ground and orbital detectors are given by relations, covering the lower 
limit, the upper limit and the tolerance of all independent test variables. 
It appears that crossed-beam missions could be accomplished with a single 
satellite, which carries two modified astronomical telescopes (ATM), which 
point toward the earth instead of away from the earth. 
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FIG. 2. REMOTE SENSING OF LOCAL CHANGES IN RADIATIVE POWER 
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FIG. 5. EVALUATION OF SPACE CORRELATIONS 
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CROSSED BEAM CORRELATION IN A SUBSONIC JET 
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OF ATMOSPHERIC RADIATION BACKGROUND 
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o o 

LU O 

o. o 

co 


o 

o 

m 


65 








LOG f MODULATION FREQUENCY (cps) 



VARIABLE DENSITY 



67 


TRACER AND AIR DENSITY DISTRIBUTION WITH HEIGHT 
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FIG. 19. GROUP VELOCITY MEASUREMENT 
WITH SCANNING BEAMS 
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APPENDIX A 

Approximation of Point Correlations in Locally Isotropic Fields 


The three dimensional wave number components of the scalar k 
field [8] can be expressed by the experimentally accessible two-beam 
product mean value R by using equation (17). 
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Sk is also indirectly related to the desired two-point product mean value 
by an inverse Fourier transformation 
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For isotropic turbulence, the wave number components are constant 
inside a spherical shell in the wave number space ^ = (^ x , fe) where 

the volume of the shell is 4 x d\|r x and the value of Sr inside the shell 
is SkOh., 0, 0, t) . Thus, the inverse Fourier transform may be expressed 
in the form 
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Substituting equation (A2) into (A3) and inverting the order of integra- 
tion gives the desired relation between the two point product mean values 
R k and the two beam product mean value R. 
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This equation can be simplified by noting that the integral can be 
expressed by the derivative of another double integral. 
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the inner integral of which is proportional to the Dirac Delta Function 
[25]. 
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Substituting (A6) into (A5) gives the final and simple result, which was 
quoted in section III, 2. 
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